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Abstract 

Tendency to bifurcate can often be utilized to improve performance characteristics of amplifiers 
or even to build detectors. Bloch oscillating transistor is such a device. Here we show that bistable 
behaviour can be approached by tuning the base current and that the critical value depends on 
the Josephson coupling energy Ej of the device. We demonstrate record-large current gains for 
device operation near the bifurcation point at small Ej. From our results for the current gains at 
various Ej, we determine the bifurcation threshold on the Ej - base current plane. The bifurcation 
threshold curve can be understood using the interplay of inter- and intra-band tunneling events. 
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In small Josephson junctions, charge and phase reveal their conjugate nature in sev- 
eral macroscopic phenomena. The quantum nature of the phase variable was shown in 
macroscopic tunneling experiments while its conjugate relationship to the charge has 



been shown in many consequent studies [2]. One of the consequences of the charge-phase 
conjugate relationship is the Coulomb blockade of Cooper pairs which arises in ultra small 
Josephson junctions having a capacitance {(7) in the femtoFarad range 0, 4]. Charging en- 
ergy and the Josephson coupling energy Ej{ip) = —Ej cosip are the competing energy scales 
associated with these two variables. Accordingly, the Hamiltonian for the small Josephson 
junction contains a periodic potential, and hence, Bloch states with band structure appear. 
These bands are analogous to the conduction electron energy states in solid state physics 

The Bloch Oscillating transistor is a three-terminal mesoscopic device which is based on 
the dynamics of the Bloch bands in a voltage biased Josephson junction (JJ) in a resistive 
environment [3, Q] . The operation is due to an interplay of coherent Josephson phenomena 
and Coulombic blockade of charge transport which is controlled by single electron tunneling 
events. The device can be viewed as a charge converter of single electrons, induced from 
the base electrode, in to a sequence of N sequential Cooper pair tunneling events, i.e., 
Bloch oscillations on the emitter terminal with a Josephson junction. The current gain is 
ideally given by /3 = 2N + 1. The number of Bloch oscillations is limited by inter-band 
transitions caused by Landau-Zener (LZ) tunneling which depends exponentially on the 
band gap between the ground and excited states of the Josephson junction. This simple 
picture has been found to correspond quite well to the measured current gain [Oj. 

Incoherent tunneling of Cooper pairs and electrons, however, complicates the basic BOT 
operation. The interaction of tunneling electrons or Cooper pairs with the electromagnetic 
environment, has been demonstrated to be strong in small tunnel junctions, both in the 



normal and superconducting states 
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Inelastic effects may, for example, limit the 
lifetime of the Coulomb blockaded state and, consequently, bias-induced changes in the 
inelastic tunneling rates can cause large modifications in the operating point, and thereby 
contribute to the current gain of the BOT. These effects, in fact, are the foundation for 
bifurcation in the BOT operation because they allow the existence of two steady states at 
a fixed base current Ib- The existence of a bifurcation point is important as, with proper 
design, the vicinity of such a point can be employed to improve the characteristics of the 



BOT. In this paper, we investigate experimentally the bifurcation threshold in the BOT, 
and demonstrate record large current gains for small- £'j-device operation near the threshold. 
From our results for the current gains at various Ej, we determine the bifurcation threshold 
curve on the Is - Ej plane. The measured transition curve can be qualitatively explained 
using a simple analytic approach, in which intra-band transitions are taken in to account 
phenomenologically, together with the transition rates due to inelastic tunneling. 

This paper is organized as follows. In Sec. I, we first outline the basic principles for 
understanding the electron tunneling dynamics in a Bloch oscillating transistor. We will 
concentrate to the dynamics near the bifurcation point at which the current gain of the 



device diverges. Our ana^ 



as done in Refs. Il2| . Sec. II will describe sample manufacture and experimental 

measurement techniques. Experimental results are presented in Sec. III. We will present 
data on the current gain at various values of Josephson energy, and construct a curve for 
bifurcation threshold on Ej vs base current plane. The relation of current gain with the 
distance from the bifurcation point is also studied in detail. In Sec. IV, we discuss our 
results in the light of analytical and numerical calculations. 



ytic model is verified using numerics with a similar approach 



I. THEORY 

A. Band model of mesoscopic Josephson junctions 

In mesoscopic tunnel junctions, the discreteness of charge starts to play a role via the 
Coulomb energy E^ = where C is the capacitance of the junction and Q is the charge on 
the capacitor plates. In quantum theory, charge is described by the operator Q = —i2e-^, 
where (p denotes the phase difference of the order parameter fields across the junction. This 
operator is canonically conjugate to 0, i.e., [Q,0] = i2e. Hence, there is a Heisenberg 
uncertainty relation, AQAip ~ 2e, which implies that the charge and the phase of the 
superconducting junction cannot be defined simultaneously. This leads to delocalization 
of the phase and to Coulomb blockade of the supercurrent, as experimentally shown by 
Haviland et al. [3] in the case when Josephson energy is on the order of the single-electron 
Coulomb energy, i.e., Ej/Eq ~ 1. The same conclusion of delocalization applies even for 



large values of the ratio Ej/Eq 113| 
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Using the differential operator due to the commutation relation, we can immediately 
write the quantum mechanical Hamiltonian as 

When Ec ^ Ej, charge is a good quantum number, which leads to Coulomb blockade of 
Cooper pairs and a complete delocalization of the phase. Equation [T] then takes the form 
of the Mathieu equation with the well-known solutions of the form = e*'^'^/^'^u„(v?), 

where Unif) is a 27r-periodic function and the wave functions are indexed according to band 
number n and quasicharge q. Verification of the existence of the energy bands has been 
carried out by different methods 14- IgI. 



Voltage over the junction is given by = ^ which changes along the energy band when 
quasicharge is varied. Thus, to have current flowing in the junction, the bias voltage Vc 
(on the collector, cf. Fig. [T]) has to be larger than the maximum Coulomb blockade voltage 
of the lowest band Eq. Vc > ^^|max- If the current through the junction is low enough, 
dq/dt ^ e6Ei/H, where SEi is the gap between the first and second band, the quasicharge 
q is increased adiabatically and the system stays in the ground band. The junction is then 
in the regime of Bloch oscillations; the voltage over the junctions oscillates and Cooper 
pairs are tunneling at the borders of the Brillouin zone, i.e., here at g = ±e. Consequently, 
the current through the junction is coherent and the voltage and charge over the junction 
oscillate with the Bloch oscillation frequency 

fB = //2e, (2) 

where / is the current through the Josephson junction. 

If the current / is not adiabatically small, we can have Zener tunneling between adja- 
cent energy bands. The tunneling is vertical, i.e., the quasicharge does not change. The 
probability of Zener tunneling between bands n — 1 and n when Ec ^ Ej is given by 

-"-^exp(-^M^e.p(-^), (3) 



where 6En = En — En-i and Iz is the Zener breakdown current [17H20|. Provided that 
Vc < ^^Imax for the excited state the junction will become Coulomb blockaded on the 
band Ei after a Zener tunneling event, and no current will flow through it any more. The 
role of the third terminal is to relax the Josephson junction back to the ground state where 
a new sequence of Bloch oscillations can be started. 



B. Incoherent tunneling processes 



The external environment gives rise to current fluctuations that couple linearly to the 
phase variable. These can cause both up- and downwards transitions. The amplitude of 
the fluctuations is given by the size of the impedance: the larger the impedance the smaller 
are the current fluctuations and the transition rates. As we will see later on, the successful 
operation of the BOT requires one to control both the upwards and downwards transition 
rates. When modeling the BOT analytically, we will make use of the Zener transition rates 
and transitions due to charge fluctuations, both derived in Ref. 3|. 

The electromagnetic environment around tunnel junctions affects the tunneling process 
by allowing exchange of energy between the two systems Q, [2l[j2^ . The influence of the 
external circuit can be taken into account perturbatively, for example, using the so called 
P(-E')-theory 23|. A perturbative treatment of the Josephson coupling term gives rise to 



a result for incoherent Cooper pair tunneling [21|, |23(] where the tunneling electron rate is 
directly proportional to the probability of energy exchange with the external environment 
governed by the P{E) function. Taking both positive and negative energy exchange into 
account, tunneling both inward and outward direction leads to the total current 

I{V) 



TTP 

{P{2eV) - P{-2eV)) . 



n 



(4) 



The function P{E) can be written as 



PiE) 



2'Kn 



dt exp 



J{t) + -Et 
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(5) 



which is the Fourier transform of the exponential of the phase-phase correlation function 



j{t) = Mt)-^m^m. 



(6) 



The phase-phase correlation function is determined by the fluctuations caused by the envi- 
ronment and it can be related to the environmental impedance via the fluctuation-dissipation 
theorem. 

For a high-resistance environment, the P{E) function is strongly peaked at energies 
around Ec-, and it may be approximated by a Gaussian function 



P{E) 



y/AnEcksT 



exp 



{E-E^ 
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^EcksT 



(7) 



where the width is governed by thermal fluctuations in the resistance R. Consequently, 
the subgap /F-curve displays a rather well-defined peak centered around V = 2Ec/e due 
to the 2e charge of Cooper pairs. This characteristic feature of the IV curve provides a 
straightforward way to determine Ej of the investigated devices of small Ej. 

The actual downward and upward transition rates Tini{Vc) and r-|-(Vc') as a function of 
the collector voltage were calculated by Zaikin and Golubev [24]. The Zener tunneling rate 
in a resistive environment, and with the assumption Ej <^ Ec, is given by 



rt = ^exp 



1 + 



and the down relaxation rate due to charge fiuctuations is given by 

f (v-l?) 



vz 



exp 



-.^27r ((5gVe2) 
where v = CVc/e, r = RqC , = ksCT, and 



2 ((5g7e2) 
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The voltage vz is related to the so called Zener break down current by Iz = evz / (4r) 



C. BOT modeling near the onset of the bistability 

Our present model generalizes the previous analytic BOT theories [s], by including 
the effect of intra-band transitions. The circuit schematics for the basic BOT modeling 
is depicted in Fig. [1] below. The basic circuit elements are the Josephson junction, or 
superconducting quantum interferometric device (SQUID) geometry, at the emitter, with 
a total normal state tunnel resistance of Rjj, the single tunnel junction at the base with 
the normal state resistance R^, and the collector resistance Rc- The BOT base is current 
biased via a large resistor Rb at room temperature, but a large line capacitance Cb results 
in an effective voltage bias. 

As required by the P(i?)-theory, our basic modeling is valid provided EjP{2eV) <^ 1. 
The intrinsic relaxation is detrimental for BOT operation and, thus, the fluctuations should 
be kept low by requiring that Rc ^ Rq = h/Ae^. In practice, we need Rc ^ 100-Rq to be 
close to the presumed idealized operation. Experimentally, this is quite hard to realize (see 
Sec. Ill) 
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Numerical analysis is needed to calculate properly the characteristics of the BOT devices 
near the onset of bistability. However, by introducing a phenomenological variable that 
describes the average number of the tunneling events {Ng) before a downward transition 
is triggered by the base electrons [12], we may derive a rather simple description for the 
operation of the BOT. A value of (A'^g) ^ 1 is facilitated by intra-band transitions that 
basically maintain the bias current of the operating point. Changes in the ratio of the bias 
current and the triggering current can lead to significant changes in the characteristics of 
the BOT. 

Like in the earlier analytic descriptions, the BOT emitter current can be thought of as the 
result of being in either of the two states; the Bloch oscillation state with a time-averaged 
constant current and the blockaded state with zero current, 

Vc/Rc, = 1/rt ^^^^ 

0, t^ = 1/{T,^ + Tb/{N,)). 

The amount of time the system spends in each state is given by the Zener tunneling rate, 
r^, the intrinsic relaxation Fin^, and the quasiparticle tunneling rate F^; only every (A''e)th 
of the injected base electrons is able to make a downward transition. The base current, 
however, flows during the opposite times: 

0, rt = l/Ft 

Ib = { (12) 

eFs, r; = l/(Fi„; + FB/(Ar,)). 

From these equations we can simply derive the average emitter and base currents 

Ce) = (13) 

Rc + n 

{Is) = (14) 

^t + n 



where we have defined 



(K) = — . ■ (15) 



By combining these two equations, we may write 



Now, when calculating the current gain Pe 
of (/_b). Thus, we obtain 

Vc Vc d{K) 

Rce{N':) Rce{N'fd{lB) 

which can equivalently be written as 

Vc {Ib) Tt(r| + Tj,) 1 



, {N'^ has to be considered function 



(17) 



with 



Pi 



Ph 



Rc e2 (AT^)^ 



(18) 



(19) 



When becomes equal to one, the gain diverges, which marks the threshold for bifurcation. 
In the regime > 1? two stable solutions are available and the operation becomes hysteretic 
as observed both experimentally and numerically. Hence, we may consider /3h as a parameter 
controlling the proximity of the bifurcation threshold. 

For ^1, we obtain a hnear dependence between and {Ib) as given by 



R, 



c 



^C7n(rt + rJ {Ib)' 

Rc + n 



{Ib) - 

n 



{Ib) 



i-B-H 



(20a) 



(20b) 



Vc n 

where e < is a phenomenological parameter to account for the variation of d (N^) /d (r^) 
under various biasing conditions (see Appendix [A]) . The latter term in the parenthesis of 
Eq. I20al specifies the threshold current Ib~h for the bifurcated, hysteretic threshold. By 
substituting (Ib) from Eq. [T3]to the prefactor of Eq. I20at {N^)'^ and {Ib)'' term cancel each 
other leaving the prefactor with {Rc/Vc){t-^ + ti)/t^. The detailed derivation of the analytic 
formulation is outlined in Appendix [Bl 

Using a simple approximation for the variation of {N'^) with r^, we may derive an analytic 
formula for the bifurcation threshold on the Ej vs (Ib) plane (see Appendix |A]) . The Ej 
dependence of Ib-h comes mainly from Eq. [H] which leads to the analytic form given by, 

Ib-h ^ Tset + exp{-KEj) 

where the first term in the numerator, Tgei^ is the upward transition rate due to single 
electron tunneling, whereas the second term arises due to LZ tunneling. The parameter k 
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involves all the other parameters inside the exponent of Eq. [SI This functional dependence 
between Ib-h and Ej in Eq. HT] is also in good agreement with the results of our numerical 
simulations. 

The BOT behaviour described here is referred to as 'normal' operation. In this configu- 
ration, the junction is initially in the upper band and a quasiparticle tunneling due to base 
current will bring the junction to the lowest band where it performs Bloch oscillations. This 
coherent oscillation will be inhibited by Zener tunneling and the system jumps back to the 
upper state and the whole process is repeated again. If the sign of Vc (and consequently 
Ie) is reversed the base current will induce transitions to the upper band, an operational 
mode that we call 'inverted' operation. Since the 'normal' operation is conceptually clearer 
we have concentrated our studies in this mode of BOT. 

II. FABRICATION AND MEASUREMENT 

The BOT samples employed in this work were fabricated using a 20-nm-thick Ge mask on 
top of LOR 3B resist. Patterning of the Ge layer was performed using conventional e-beam 
lithography at 20 keV. After patterning, the PMMA layer was developed in MIBK:IPA (1:3) 
solution and subjected to a plasma etch with CHF^ plasma. Finally, the LOR under the 
germanium was etched in oxygen plasma up to the desired extent of undercut. 

Shadow angle evaporation at four different angles was employed to generate the structures 
consisting of three metals. Originally, the BOT was envisioned to have a NIN junction as 
the base junction, but the technique of manufacturing both SIS and NIN junctions on the 
same sample is exceedingly difficult and, therefore, we opted to have a NIS base junction 
instead. The SIS junction is formed of two Josephson junctions in the SQUID geometry; 
this facilitates tuning of the Josephson energy by magnetic fiux. The process order in the 
evaporation sequence was (I) Chromium, (II) Aluminium, (III) oxidization, (IV) Aluminium, 
and (V) Copper. NMP or PC remover was used for lift-off. Oxidation was done in Ar:02 
(6:1) mixture at 80 mTorr for 1 min. 

A typical sample used in the present study is displayed in Fig. [H The area of the SIS 
junctions is 100 x 150 nm^ each (equal areas within 10%). The NIS junction on the base 
has an area 70 x 100 nm^, roughly half of the SQUID junctions. 

The measurements were done on a plastic dilution refrigerator (PDR-50) from Nanoway 

9 





FIG. 1. Scanning electron micrograph of the sample (left frame) and a schematic view of the device 
(right frame). In both pictures, base, emitter and collector are marked by B, E and C, respectively. 
Positive directions for the currents are indicated by the arrows. The sample parameters are given 
in Table m Qi{t) is the island charge tracked in the numerical simulations. 



BOT # Rn Rjj Rc Ej Ef^ Ec A 

1 53 27 550 17 2.7 40 150 

2 75 21 305 25 3.3 60 165 



TABLE I. BOT parameters for the measured sample. Rm and Rjj are the normal state resistances 
of the NIS and JJ tunnel junctions in the SQUID-loop geometry, respectively. Resistances are given 
in units of kO and energies in /ieV. 



Ltd. The base temperature of the refrigerator was 50 mK. The filtering in the PDR consisted 
of 70 cm long Thermocoax cables on the sample holder and 1 kOhm serial resistors at 1.5 K. 
In addition, micro- wave filters from Mini-circuits (BLP 1.9) were used at top of the cryostat. 
The measurement set-up in this work was similar to that described in Ref. |26|. The 



BOT base was DC current biased by a resistor Rb=1 — 10 GQ, which was located at room 
temperature. Voltages were measured with low noise LI-75A voltage preamplifiers while 
currents were monitored using DL1211 low noise current amplifiers. 

The resistance values of the three circuit branches were determined at 4.2 K. Since 
there was a weak temperature dependence in Rc, we determined the actual value from 
1/VV asymptote 27 1 of the /F curves measured at low Ej. The maximum Josephson 
energy Ej was calculated using the Ambegaokar-Baratoff relation which yielded Ej = 
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17 /ieV. The flux-modified Josephson energy was obtained from tlie formula = 
-E'jYCos^(7r$/$o) + (i^ sin (7r$/$o); where d = ^^"^ denotes the asymmetry in Joseph- 
son energies between the two SQUID loop junctions with Ej^ and Ej^^ respectively. By 
fltting Ej{(f) to the measured IV curves we found d = 0.15 and 0.13 for samples #1 and 
#2, respectively. Emitter-collector and base-emitter /F curves were employed to determine 
the effective energy gap A of the samples (see Table [T]), which is 20 — 30 /ieV smaller than 
the bulk value A = 0.18 meV. This reduction of the gap is presumably due to the inverse 
proximity effect [28] due to the chromium resistors. The spatial variation of the inverse prox- 
imity effect would also explain the larger asymmetry between the SQUID junctions than is 
expected due to the difference in their areas. 



III. EXPERIMENTAL RESULTS 



IV characteristics 



IV characteristics of sample ^1 measured at a few magnetic flux values are illustrated 
in Fig. [2J the emitter-collector current Ie is recorded as a function of at = 0. The 
data clearly shows Coulomb blockade of supercurrent at all investigated values of 

the Josephson coupling energy. The peak in the IV in the subgap region is a signature 
of the inelastic Cooper pair tunneling, commonly referred to as P(E) peak (cf. Eq. [7j). 
The weakness of the blockade in Fig. [2] is assigned to the small Coulomb energy Ec = 40 
fieV, the value of which is determined from the position of the P{E) peak, extrapolated to 
Ej = 0. At larger bias voltage, Zener tunneling to higher bands takes place, which causes 
the phase fluctuation theory to break down. Our results on Zener tunneling are similar to 
those of Kuzmin et al. who investigated a single Josephson junction in an environment of 



chromium resistor 



3- 



Fig. [3] demonstrates the effect of the base current on the IV of the BOT. The 'normal' 
and 'inverted' operation regions are deflned by sign combinations {Vc, Ie, —Ib) and (Vc, Ie, 
Ib)i respectively Both the 'normal' and 'inverted' modes of operation display a strong 
increase in the onset of the LZ tunneling current, which is seen as the movement of the 



shoulders in the IV curves up to larger currents. The down-turning s 
and lower sets of the IV curves are distinct features of LZ tunnehng 
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loulders in the upper 
while the data at 



-0.5 0.0 0.5 

V, (mV) 



FIG. 2. /F characteristics of sample ^^1 at a few values of Josephson coupling energy: Ej = 9 fieV 
(black), 7.3 /ieV (red), 6 /ieV (blue), 4.5 /ieV (magenta), and 2.8 /ieV (green), without base current 
{Ib = 0) at T = 90 mK. 

Ib = display only smeared bumps of these features. This enhancement of the LZ current 
suggests that effectively the energy gap between the ground and excited states is increased 
due to the noise induced by the current in the base junction. The shoulders move even 
further apart with growing base current, which indicates an increase in the effective energy 
gap at the Brillouin zone boundary. In general, the 'inverted' operation displays comparable 
characteristics as the normal operation, but we found that hysteretic behaviour appeared 
at smaller bias currents in the inverted regime compared with the normal operation mode; 
in some cases, these modes differed by a factor of four in the base current for bifurcation 
threshold. Nonetheless, since the normal operation is appeared to provide more clear-cut 
data, we concentrated our studies on this operating regime. 

B. Gain determination 

Fig.lSb displays a basic set of data for current gain determination in the 'normal' operating 
region. Emitter current Ie is depicted as a function of collector voltage Vc at eight values 
of base currents Ib- The regime with a large negative slope marks the active bias regime 
of the BOT amplifier. The steepest monotonic curve (the second one from left in Fig. [3]d) 
has a narrow linear regime in the center of the negative slope part, the width of which 
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FIG. 3. (a) Middle traces: magenta and black IV curves are measured without base current at 
Ej = 6.5 and 5.8 /ueV, respectively. Red curves, corresponding to Ej = 6.5 //eV are measured 
at Ib = +0.3, -1-0.34, and -1-0.38 nA (traces from right to left). Blue curves have the same bias 
conditions as the red curves but were measured at Ej = 5.8 ^eV. The red curves are offset by 
(-F0.22 mV, -hO.42 nA) for clarity, Hke the blue curves by (-^0.22 mV, -0.42 nA). 
(b) The normal operation region of the BOT at Ej = 7.1 /ieV with increasing Ib- Negative 
slope is the Landau-Zener tunneling regime, increases with Ib and eventually the slope diverges: 
Ib = -h0.06, 0.065, 0.07, 0.075, 0.08, 0.085, 0.095, and 0.105 nA (from right to left). Filled 
(open) circle traces are of Ie when Vc is swept from left (right) to right (left). The measurement 
temperature was at T ~ 90 mK. 



amounts to about 2 pA in Ib- This corresponds to the maximum dynamic range in Ib over 
which the BOT has substantial current gain at this bias point. Roughly, a change in the 
base current by A/b = 2 pA corresponds to 50 pA in Ie-, and the current gain becomes 
Pe = 25. Eventually, the slope of the IV diverges with increasing Ib-, after which the IV 
characteristics become hysteretic as seen at the largest value of Ib = 0.105 nA in Fig. [Jb- 
Clearly, in the BOT operation near the divergence point, the dynamic regime is inversely 
proportional to the current gain. 

We have checked that hysteresis does not depend on the value of the current bias resistor 
in the range 10*^ — 10^° Vt. Moreover, we have performed simultaneous transconductance 
Qm = and current gain measurements to determine the input impedance of the BOT 
■^^"^ ~ ~ ^ ^if ~ g^- ^^^^ input impedance diverges at the same point 
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0.04 0.06 0.08 0.1 0.12 



FIG. 4. Measurement of the current gain by tracing Ie vs Ib at Ej= 7.1 /^eV. The steepest slope 
yields the operating point with the largest current gain (3e at the corresponding collector voltage 
Vc- Traces were measured at T = 90 mK using Vc = -0.443, -0.429, -0.419, -0.410, and -0.401 
mV (traces from right to left). Different signs of Ie and Ib correspond to the regime of 'normal' 
operation. Red (purple) traces are for growing (decreasing) sweep of Ib- 

as the gain. 

According to basic BOT theories [3, Q, 25 1, the current gain is independent on the base 



current. However, the situation changes near the bifurcation point. This is because there 
can be two different kinds of base current components: one comprising of tunneling events 
causing interband transitions (the only component in the traditional BOT base current) and 
another one leading only to intraband events. Only the interband transitions lead to gain 
in the BOT while the intraband transitions are to maintain the bias current. As the base 
current grows, the ratio of these two current components may change with increasing Ib 
which leads to current dependence of the gain and, eventually, to the diverging behavior 
when approaching the bistability point. Hence, the observed strong increase in Pe with 
increasing base current is a sign of the operation near the bifurcation point where the gain 
grows according to Eq. EDI 

In addition to the analysis of data like in Fig. [Hb, we have measured the current gain 
using traces of Ie vs Ib as illustrated in Fig. HI The figure displays data at five different 
values of the Vc at Ej = 7.1 yueV. The current gain is calculated from the negative slope 
of Ie{Ib)'- Pe = "^if • The maximum of the decaying slope yields the optimum current 
gain, which we determined as an average of the up and down Ib sweeps. At large gains, 
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there was often a rather large difference (~ factor of 2) between the gains of up and down 
sweeps. In such cases, we considered the operation of the BOT bistable at this bias point and 
disregarded the larger gain value. We performed Ie{Ib) measurements at different values of 
Ej and picked out a single gain value I3e(Ib) = ~X7^ from each of the data sets. When 

max 

Pe > 50, the gain determinations became problematic because of 1// noise and creep in 
the measurement, which gradually took the device out of the linear regime during the gain 
determination. 

Fig. O depicts data on the inverse of Pe vs Ib which were obtained from the analysis 
of Ie{Ib) scans performed in the range with Ej = 4.6 — 10.5 /ieV. Plotting (3^^ makes the 
analysis of the diverging gain regime simpler and allows us to examine the vicinity of the 
bifurcation point where we are supposed to have (3^^{Ib) 0. Experimentally, the problem 
in this region arises because the dynamic range becomes zero and measurements without 
noise induced smearing become impossible. Nevertheless, the data in Fig. O display how the 
critical regime is approached at (3^^ > 0.02 which corresponds to our highest reliable gain 
values. All of the data at small values of (3^^ are seen to show a nearly linear dependence on 
Ib, especially at large values of Josephson energies. At our smallest value of Ej = 2.7 /xeV 
(sample #1), we could not reach the bistable regime at all. The theoretical dependence for 
f^E^i^s), illustrated by red curves in Fig. were obtained by fitting Eq. I20al to the data 
just near the divergence point, as required by its regime of validity. 



C. Bifurcation threshold 

The experimentally determined values of Ib-h for the bifurcation point are plotted in 
Fig. Elon the Ib-h — Ej plane. The plot was generated from the fits in Fig. \5\hj selecting the 
points of (3^^{Ib-h) = 0. Fig. [H] indicates that the onset of bistability is nearly independent 
of base current at large values of Ej, while a steep increase in Ib-h is observed below Ej = 6 
fieV. The observed behavior is quite well reproduced by our phenomenological formula Eq. 
[2T] which is depicted by the solid curve. 

For sample 7^2, we found a similar threshold curve which indicates that the bifurcation 
behaviour and its dependence on Ej is a fundamental characteristic of Bloch oscillating 
transistor. The bifurcation threshold curves for both the samples are depicted in Fig. O 
For sample #2, the bifurcation threshold current is higher than that for sample #1. From 
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FIG. 5. Inverse gain as a function of bias current 1b- Each data point was obtained from 

a 1e vs 1b sweep illustrated in Fig. HI The sohd curves were obtained using Eq. [2D] fitted to the 
highest 1b quartile fraction of the data sets (lst-8th lowest values). 



the fitted curves we found that k (see Eq. I^Tl) for sample #1 is higher than for sample 
which comes from the fact that k contains term Rc which is higher in sample #1 than in 
sample #2. The smaller base current needed for bifurcation for sample #1 than sample #2 
is accounted from the fact that fin for sample #1 is higher than for sample #2. Absence of 
bifurcation was observed in both sample #1 and #2 at their respective lowest Ej values. 



The rate 



at which l//?^; reaches the bifurcation point depends on Ej as seen from 



Fig. 7(a) Initially, the slope increases rapidly with Ej upto 6 /ieV, while between 6 - 10.5 
/xeV the slope appears to be saturated. In this region, the bifurcation threshold current is 
almost indepedent of Ej (cf. Fig. E]). The slope variation over the whole range of Ej amounts 
to a factor of 5.5. Hence, the variation of Ej does not necessarily change fiE strongly, which 
is a desirable property concerning 1// noise due to critical current fluctuations. 
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FIG. 6. Bifurcation threshold on the Ej vs plane. Red(filled) and blue(open) circles denote 
the Ib-h values for the sample ^^1 and ^^2, respectively. Solid curves display the theoretical 
dependence from Eq. [2TJ 
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FIG. 7. (a) 



VS Ej plotted near the bifurcation threshold. Each point at different Ej was 



derived from fits in Fig. [5l (b) deduced from experimentally determined p^r^^, is plotted 

vs Ej. Solid curves are seen to display the similar character but the theoretical estimate falls short 



by 70% from the change in Fig. 7(a 



Theoretically, the rate ( ) is hard to evaluate from Eq. [2D1 The significant pref actor 

^ AO — 1 

of Eq. I20bl contains two terms : Vc/Rc and (r^ + r-t-)/r|. We have analyzed how varies 
theoretically with Ej for sample #1 by estimating the factor By using Eq. [HI 

we can relate -^"^p- to the experimentally determined quantity y^j^^^ , while Vc is obtained 
from the bias voltage, which is increased by 60% over the range of Ej = 3.9 — 10.5 /xeV. In 
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Fig. 7(b) , we show the variation of -^"^^ with Ej as determined for the maximum current 
{Ie) at the sub-gap I{V) peak. Together, these opposing contributions result in a change 
by a factor of 3.3 in , which falls slightly short from the observed factor of 5.5 in Fig. 



7(a) Hence, we can conclude that our simple model explains the rate of approach towards 



the bifurcation threshold with fair extent. 

We have also tried to determine the ratio of inter- and intra-band transitions which is 
governed by (A^e)- According to Eqs. [T3] and dS] there is the following relation between the 
base current and (Ne) : (/b) = ^^^^ T-t+n " ^^^^ approximation, {Ib) oc F^, we can 



conclude that {N^) oc Is-Hy ^^^'^ ■ From Fig. 7(b) the increase in ^'^^^^ is ~ 2 whereas 
from Fig. [H]the decrease in Ib^h ^oi sample #1 is ~ 4. Thus, we can conclude that (A^e) 
goes down with increasing Ej. But unfortunately, we cannot determine the exact number 
of (A"e) from this analytical formulation. 

In our numerical analysis, we have considered the circuit model used by Hassel et. al js], 
and modified it for current bias configuration. Here we have calculated the island charge 
as a function of time by taking into account three contributions: charge relaxation through 
Rc together with the tunnel current through the emitter and base junctions, respectively. 
The tunnel currents through these junctions are calculated using time dependent P{E) - 
theory Q]. In the simulation, P{E) is calculated numerically by considering only the real 
part of the environmental impedence. The simulation-run-time was chosen longer than the 
time constant due to Rb and the capacitance from base to ground so that the steady state 
was reached properly. Moreover, we monitored the tunneling events on the island with time, 
which clearly revealed the Bloch oscillating state and its transition to the higher band. By 
counting the number of tunneling events when the system undergoes a change from higher 
band to lower band we could calculate (Ng) from the simulation. 



IV. DISCUSSION AND CONCLUSIONS 



According to the basic theory of the BOT operation, the gain should depend exponentially 



on Ej/Ec via the tunneling rates F-i^ and Tmi [25|| when Ej <^ Ec- In this small- E'j limit, 
the energy gap between the first two bands is small, which facilitates the use of the up- and 
down transition rates (Eqs. [S] and ^ from perturbation theory. In our experiments we are 
well in this limit, which has not been the case in many of the previous measurements, for 
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example, in Ref. Q the maximum gain of Pe = 35 was achieved for Ej/Ec = 3.4. A current 
gain of Pe = 25 was reported for low Ej/Ec = 0.3 in Ref. y. In our present paper, we 
have observed a large current gain of ~ 50 even at Ej/Eq = 0.1. The estimates from the 



basic theory 25| amount to /3e = 4.8 — 7 for Ej = 5 — 11.8 /ieV, well below the measured 
values. Moreover, we did not observe any variation of the maximum gain with Ej, which, 
together with magnitude of Pe, is consistent with the operation near the bifurcation point 
where the main gain mechanism has a different origin than in the regular BOT operation. 

In the operating regime near the bifurcation point, the base current is a combination 
of a working point current, not inducing inter-band transitions, and a significantly smaller 
part that leads to transitions, the ratio of these two currents being given by the parameter 
(Nf.). In our phenomenological modeling with large number of intra-band transitions, the 
current gain is simply related to {N'J) and the upward tunneling rate (see Eqs. [15] and 
[T6i) . Hence, a large current gain is expected when approaching a regime where there are two 
stable solutions for the base current with different values for (A^e)- When (A^^e) is large, then 
almost all of the current is used to just keep the operating point. In our simulation we find 
a factor of 15 change in (Ag) over the measured range but unfortunately, our analysis is not 
able to yield absolute numbers for (A'g) from our measured data. 

When comparing our findings with the numerical work of Hassel and Seppa , we 

find a weaker overall dependence of the device performance on the sample parameters and 
biasing parameters than was found in the simulations. The weaker overall dependence may, 
of course, be valid only for the regime of the sample parameters/device configurations that 
were investigated in the present work. Nevertheless, the weaker parameter dependence is 
an important factor that contributes to the success of the simple phenomenological mod- 
eling that we have employed. The weaker overall changes may also indicate that there 
is external noise present in the measurements and our results should be compared with 
simulations performed at a higher effective temperature. Furthermore, as an example of 
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the differences, let us point out that if we take the hysteresis parameter from Ref 
Ph = 0.02 exp (^ ^'^^^f (^) ) ' ^^^^ sample #1 should be bifurcated at all 

base currents {Ph ~ 2.4 — 11.5). We think that the absence of bifurcation at = with 
Ph well above 1 indicates the necessity to add a capacitance in parallel to Rc into the sim- 
ulations, which would take into account the parasitic capacitance component on the sample 
chip. This parasitic capacitance will influence the Coulomb blockade at large frequencies, 
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at which it will reduce the real part of the impedance seen by the Josephson junction. 

In conclusion, we have investigated the dynamics and modeling of the BOT when ap- 
proaching a bifurcation point governed by intricate interband transition dynamics. Our 
results present the first experimental analysis on the behavior of the BOT in the regime 
where its behavior is fully governed by switching dynamics with the rates imposed by the 
biasing conditions. We have reached record-large current gains even though the device was 
operated just at small Josephson coupling energies Ej = 2.7 — 10.5 /ieV. We have mapped 
a cross-over transition diagram on the Ej vs Is plane and compared its shape to analytic 
modeling, where the intraband transitions are included in terms of a phenomenological pa- 
rameter (Ne). The same modeling was also successfully applied to describe how the current 
gain diverges as a function of the base current Ib- Our findings are consistent with the gain 
divergence as 1/{1-(3h) where the hysteresis parameter (3h is only weakly dependent of Ej. 
The weak dependence makes this regime attractive for application where large current gain 
is needed at low frequencies. 
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Appendix A: Analytical derivation of threshold curve 

It is difficult to obtain an analytic expression for the derivative d (N^) /dr^ and hence, 
we had to be satisfied with crude approximations. Eq. ^\ specifies the relation between 
I3h and the derivative d {N'^) /dr^ as follows: (3h = jf^^^~'- The partial derivative of {N^) 
with respect to can be approximated as 

^^<^ + .. (Al) 
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where e < is a phenomenological correction term. In order to determine the variation of 
{N'^) with Tj, exphcitly, we follow an interpolative approach using 

1/2 



(A2) 



which agrees with the limits; when is short {N'^) approaches 1 and (A^^^) ~ F^r^ when 
is long. Hence, 



which leads to, 



9n n 



dn 



1 



Using Eq. lAll we can write for the hyteresis parameter, 

e d{Ni) 



/3. 



H 



iK) 



(Ib) 



+ e 



Using the expression of {Ib) from Eq. [HI 

{Ib) = 



e{K) 



(1 + ?)' 



and substituting it with -^^^^ into Eq. IA6| we get the following equation 



(A3) 



(A4) 



(A5) 
(A6) 



^B-H U 

~ — e— 



(AT) 



e 

at Ph = 1- By inserting e from Eq. IA4[ we obtain an analytic expression for the bifurcation 
curve, 



1 



1 



(A8) 



^1 + {Tsny 

The upward transition rate I/t-^, depends on both LZ tunneling and single electron tun- 
neling (rse-|-) . Single electron events induced by the base current were found to be important 
in the simulated time traces of the island charge at high Ej values. Hence, 

1 



F-f + F^et 



(A9) 
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According to LZ tunneling (F-i^) (cf. Eq. [H]) 



(AlO) 



where the parameters inside the exponent are absorbed in k. Thus, we arrive at an expo- 
nential dependence of Ib-h with E"^: 

Ib-h r^et + exp(-K^3) 



oc 



1 + rl/ Ej 



which is Eq. |2I] in the main text. The effect of single electron tunneling reflects on the 
bifurcation threshold curve through the saturation of Ib-h{^ 0) at higher Ej values. 

We have used Eq. |5T] to fit the bifurcation threshold diagram. We find good agreement 
with both experimental and simulated data (cf. Figs. M and [H]). Though our simulated 
threshold values deviate from our experimental Ejs, we find similar functional tendency in 
the curves. Both in the simulation and experiment we observe that the bifurcation takes 
place earlier in 'inverted' regime than in the 'normal' operation. In the inset of Fig. [HI 
we display a calculated (3^^ vs Ib plot at Ej=10 /leV. In the simulation, we also found a 
minimum Ej below which there is no hysteresis. Hence, we can conclude that our simulation 
quite well explains the experimental findings. 
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FIG. 8. Bifurcation threshold on the Ej vs Is plane obtained from the simulation. The soUd curve 
is the analytic dependence from Eq. \TI\ The inset shows the dependence of (3^^ on Ib obtained 
from the numerical simulation at Ej=W /xeV; the fitted line indicates Ib-h = 0.16 nA. 
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Appendix B: Is vs (Se 



According to Eq. 



_ Rc {N'f ^ o . 



(Bl) 



By inserting Ph from Eq. 
point: 



we reach the approximate form of Pe near the hysteresis 



Vc {Ib) n(n + n) 



Rc e2 

Rc it + n 



(Ib) 



-1 



Vc n 

where Ib-h denotes the bifurcation threshold current. The above formulation is valid only 
in the vicinity of the divergence point, where the dominant change in can be viewed as 
linear in {Ib)- 
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